Universal Spectra of Coherent Atoms in a Recurrent Random Walk 



R. Pugatch,^ O. Firstenberg,^ M. Shuker,^ and N. Davidson^ 

^Department of Physics of Complex Systems, Weizmann Institute of Science, Rehovot 76100, Israel 
^Department of Physics, Technion- Israel Institute of Technology, Haifa 32000, Israel 

The probability of a random walker to return to its starting point in dimensions one and two is 
unity, a theorem first proven by Polya [1] . The recurrence probability — the probability to be found 
at the origin at a time t, is a power law with a critical exponent —d/2 in dimensions d — 1,2. We 
report an experiment that directly measures the Laplace transform of the recurrence probability in 
one dimension using Electromagnet ically Induced Transparency (EIT) of coherent atoms diffusing 
in a vapor-cell filled with buffer gas. We find a regime where the limiting form of the complex 
EIT spectrum is universal and only depends on the effective dimensionality in which the random 
recurrence takes place. In an effective one-dimensional diffusion setting, the measured spectrum 
exhibits power law dependence over two decades in the frequency domain with a critical exponent 
of —0.56 ± 0.01. Possible extensions to more elaborate diffusion schemes are briefly discussed. 



A remarkable theorem due to Polya [1] asserts that the 
probability of a random walker to return to the origin at 
dimensions d = 1, 2 equals one. Such walks are called re- 
current walks. The Polya theorem has far reaching con- 
sequences in physics. For example, it is related to the 
well-known result that, in the non-interacting zero-spin 
Anderson model, disorder is strongly localizing at dimen- 
sions less than three [l, Strikingly, it can be shown 
that the absence of spontaneous symmetry-breaking 
phase-transition m. d < 3, the so called Hohenberg- 
Mermin- Wagner theorem, is also a consequence of Polya's 
theorem [4] . Given the probability to be at the origin at a 
given time P(r = 0, t), it can be shown that the walk is 
recurrent if and only if P{0^t)dt diverges [5[. Closely 
related to P(0, t) are the distribution of First Passage 
Times, FPT(r, t), which is the probability to reach the 
point r for the first time at t, and the distribution of First 
Return Times, FRT(t) = FPT(r = 0,t). These distribu- 
tions are of great interest in many fields of research, from 
reaction-diffusion processes, the onset of firing bursts of 
neurons, porous media structure analysis, transport in 
disordered and scale invariant media, and the spreading 
of diseases [El, In particular, FRT distribution de- 
termines trap performance in velocity-selective coherent 
population-trapping schemes f^. 

Here, we report an experiment that directly measures 
the Laplace transform of the one-dimensional recurrence 
probability of a Rb atom diffusing in N2 buffer gas. This 
is achieved by measuring the spectrum of Electromagnet- 
ically Induced Transparency (EIT) of a narrow light sheet 
(Fig. [1]). The complex spectrum thus obtained has a uni- 
versal power-law form, which is the Laplace transform of 
the recurrence probability, P(0, s) cx P = 0.5, where 
s = zA + ro/2 is the complex Raman-detuning, with A 
being the two-photon Raman-detuning, and Fq the EIT 
natural line-width. The measurements in the universal 
regime, where the power broadening and the beam size 
are small, depends only on d — the effective dimension- 
ality of space. Our measurement yield this power law de- 
pendence with a critical exponent of P = 0.56 ±0.01 over 



two decades in the frequency domain. In a control exper- 
iment, using a wide- area circular beam, the spectrum is 
expected to have power-law dependence of = 1, 

which corresponds to a Lorentzian spectrum. In this ex- 
periment, our measurement yield = 0.97 ± 0.01. An 
explanation of these results is given using a sum over his- 
tories approach, which can be generalized to other sys- 
tems where the atomic motion is not restricted to normal 
diffusion and can be either of the anomalous diffusion 
type or ballistic (in chaotic, mixed, or regular systems). 
An alternative derivation is also given by directly solving 
the relevant diffusion equation for the atomic coherences 
0. Our results are an extension of the pioneering work 
of Xiao et. al. [9], who recognized the importance of 
coherent returns as a line-shape narrowing mechanism, 
albeit not in the universal regime. 

The probability of a diffusing particle that starts at the 
origin at t = to be at position r at time t is given by 
P(r,t) = (47rP)t)-^/2e-^ ^^^K It follows that the recur- 
rence probability is P(r = 0, t) = (47rl)t)~^/^, which is a 
power law with an exponent —d/2. It is well known that 
the first passage time distribution is related to P(r, t) 
and P(0,t) by [5]: 



P(r, t) = Jr,o^^,o + / dt'YVT{v, t')P{0, t 
Jo 



t'). (1) 



This equation can be simply interpreted by noting that 
the probability to reach a point r at time t is the sum 
over all intermediate times t' of the probabilities to reach 
r for the first time at t' and then to return to this point 
after a time t — t'. The first term in Eq. ([T]) is a boundary 
term for t = 0. To calculate the FRT, we substitute r = 
in Eq. ([T]), take the Laplace transform from both sides, 
and recover 



YRTd{s) = 1 



Pd(0,5)' 



(2) 



where Pd{0,s) is the Laplace transform of the recur- 
rence probability in dimension d. For (i = 1, Pi (0,5) = 



2 



(ADs) which is a power law in 5, while for d = 2 
P2(0,s) = K^^^JJJVb) - ln(5/rz))/4 - 0.577 for small 
5, where is the modified Bessel function of the second 
kind, Vb — 4D/Wo is the diffusion broadening, D is the 
diffusion coefficient, and Wo is the beam diameter. 

In brief, our EIT is obtained with two phase-coherent 
laser fields that couple the two magnetically insensitive 
hyperfine ground-state levels of ^''Rb: |1) = |5^5'i/2, F = 
l,m = 0) and |2) = \^'^Si/2,F = 2, m = 0) to a mu- 
tual excited level |3) = |5^Pi/2,F' = 2, m = 1). When 
the energy difference between the two beams is equal 
to the ground level splitting, the atom is pumped to- 
wards a dark state \D) = 2~^/^(|l) - |2)) at a rate Fp, 
which is proportional to the laser intensity. An atom in 
a dark state decouples from the laser fields, and is there- 
fore transparent. An orthogonal state to \D) is the bright 
state \B) = 2-^/2(|l) + |2)), which can be excited to |3). 
An atom in a dark state that wanders out of the EIT 
beams will perform dark-to-bright oscillations at a rate 
determined by the two-photon Raman-detuning, A, in 
the rotating frame defined by the two lasers. 

An atom in the dark state that stays within the beams 
for a time t has a probability Re(e~*^^*) to be in a dark 
state upon leaving, where Sin = — + (Fq + Fp)/2. An 
atom in a dark state that leaves the EIT beams for ex- 
actly T seconds has a probability Re(e~*^) to be in a 
dark state upon returning, where s = —iA + Fo/2 [iq| . 
If the pumping rate back to dark state, Fp, upon return- 
ing to the beam, is much smaller than the dark-state 
decay rate, Fq, the total transparency is proportional to 
Re(e-^^). 

In our experiment, the EIT beams are localized to a 
small volume in space. The atoms diffuse in and out of 
the beam randomly due to the collisions with the buffer 
gas. The steady-state EIT spectrum is a sum over all pos- 
sible histories of dark-state atoms being within the beam 
for exactly t seconds, then staying outside the beam for 
exactly T seconds, then coming back for exactly f sec- 
ond, staying outside for exactly T' seconds, and so on. 
For example, the contribution to the spectrum of atoms 
that were pumped to dark state after passing through 
the beam and then lost their coherence outside the beam 
is S{s) = dtFFT{t)e-'^-^ = FPT(5in), where FPT(t) 
is the average FPT distribution through a beam of diam- 
eter Wo [11]. Summing over all histories, we obtain 



Sis) = ^FPT(si„)"+iFRT(s)" = - 



FPT(si„) 



n=0 



FPT(sin)FRT(s)- 
(3) 

Note that long trajectories are less probable due to the 
decay factor. 

In the limit when the beam size tends to infinity with 
uniform intensity, there are no dark periods and thus 
the FRTs drop from Eq. (|3]). The obtained complex 
spectrum is a power law, S = FPT(5in) = s-"^, which 
is the well-known complex Lorentzian spectrum, i.e., the 



absorptive part of is the Forentzian distribution. The 
width of the Forentzian is Fo + Fp, where Tp is recognized 
as the power broadening. 

In contrast, the universal limit is obtained when the 
beam size tends to zero Td ^ Fq, with a uniform and low 
intensity, such that F^ <C Fq. In this limit, the transit- 
time broadening is large and FPT(sin) 1. The depen- 
dence on the transit time drops from Eq. (|3]), and the 
spectrum becomes universal: 



5(s) = 5]FRT(s)" = - 



FRT(s) 



= P{0,s). (4) 



We now turn to an alternative derivation of Eq. (|4]). 
In [12] we demonstrated that the evolution of the atomic 
coherence field during storage-of-light is governed by a 
diffusion equation. In [8], we developed a comprehensive 
theory for EIT spectra of thermal atoms. The result- 
ing diffusion equation for the coherence and the complex 
spectrum that directly depends on it are given by 



-zA + ^ - DVi 



(r) 



i^2i (r) = -noFp (r) , (5) 



S{s) (X J d^r Tp (r) [2i?2i (r) /no] ; 

where Fp(r) is proportional to the intensity profile, 
R21 (r) is the steady-state ground state coherence, no 
is the density (~ 9 • 10"^^ /cc), and S{s) is the normal- 
ized complex spectrum. In the universal limit, we ne- 
glect Fp(r) from the left hand side of Eq. ([5]) and replace 
Fp(r) Tp5{r). We then obtain 

[-iA + ^ - DVi]i?2i = -noTpSir), (6) 

Taking the Fourier transform, algebraically inverting the 
equation, and taking the inverse Fourier transform, we 
find 



i?2i(r) = -noTp / d-^q 



s + Dq^' 

5(s)a / d''ri?2i(r)<5(r) =Pd(0+,s). 
Jo 



(7) 



This equation can be solved exactly in dimensions d = 1 
and d = 2. In d = 1, S{s) ex and for d = 2, 

S{s) (X Ko{^/s/^D) ^ ln(5/FD)/4 - 0.577 for small s, in 
agreement with the sum-over-histories approach [Eq.(|4])]. 
This agreement stems from the fact that one can always 
write the diffusion propagator as a path integral. Tak- 
ing the beam size to infinity while keeping the intensity 
finite, we can neglect the diffusion by setting D = in 
Eq. (5). The resulting complex spectrum is the complex 
Forentzian spectrum, S{s) cx s^^. 

Our experiment consists of a single mode Vertical Cav- 
ity Surface Emitting Faser (VCSEF) tuned to ~ 795 
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FIG. 1: Measured signal vs. the two-photon detuning A and 
the power law fit. Inset (a): zoom-in of the measured signal 
over the range =b300 Hz. Inset (b): an illustration of the one- 
dimensional beam experiment. Returning atoms contribute to 
the total intensity irrespective of the location of their reentry 
point thus effectively reducing the dimensionality to one. 
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FIG. 2: Measured signal for the Id beam configuration (tri- 
angles), and for the wide area beam (circles) vs. the com- 
plex detuning size \s\ in units of To/ 2. The lines are the- 
ory with a single fit-parameter — the critical exponent. We 
find (3 = 0.56 lb 0.01 over two decades in frequency and 
(3' = 0.97 ± 0.01 over two decades. 



nm on the ^^Rb Dl transition. By modulating the VC- 
SEL's current we created two equally powered sidebands 
at ~ ±3.4 GHz, which consist our co-propagating EIT 
modes (separated by the hyperfine splitting of ~ 6.8 
GHz). A secondary weak frequency- modulation (FM) 
at 97 kHz was added to facilitate low noise lock-in detec- 
tion. Using weak (1%) modulation depth ensured that 
the measured beat signal is directly proportional to the 
dispersive part of the EIT spectra [1^ . The beam was cir- 



cularly polarized and passed through an isotopically pure 
^^Rb vapor cell with 10 Torr of N2 buffer gas {D = 10 
cm^/sec). The cell dimensions were 2.5 cm in diameter 
and 5 cm in length, water-heated to ~ 48C, and magnet- 
ically isolated using three concentric /i— metal shields. A 
set of Helmholtz coils were used to provide a homoge- 
neous longitudinal magnetic field in the beam propaga- 
tion direction of 49 mG, thus pushing the magnetic sensi- 
tive EIT lines ±69 kHz sideways. The line- width broad- 
ening due to the inhomogeneous second order Zeeman 
shift is less than 1 Hz, much smaller than the line- width, 
which is essential for observing the universal regime. Af- 
ter exiting the cell, the EIT beam was focused on a photo- 
detector (PD). The PD signal was filtered and amplified 
using on-board electronics and was fed to a dual-channel 
lock-in amplifier (SRS model SR-850) to allow for de- 
modulation of the 97 kHz EIT beat signal. 

In our measurements, we scanned the two-photon de- 
tuning A in the range ±8 kHz in frequency steps of 2 — 4 
Hz. The dwell time per detuning was 0.5 sec. After a 
proper settling time of typically 10 msec, we sampled and 
subsequently averaged the lock-in signal and obtained the 
dispersive part of the EIT spectrum as a function of A. 
In the first experiment, we used a wide area beam with a 
diameter of 8 mm. We measured the spectrum for differ- 
ent powers and found the power broadening slope to be 
165 Hz/(W/m^) and Fq = 45 Hz. In the second experi- 
ment using a spherical and a cylindrical lens, we shaped 
the beam to become a light sheet with beam diameters of 
126 /im xl cm along the cell length [Fig. [I{b)]. Atoms 
returning to the beam contributes to the transparency ir- 
respective to the location of their reentry point [14] mak- 
ing this geometry effectively one-dimensional. The diffu- 
sion broadening, determined by the narrow dimension of 
the beam Wo = 126/im, is Fi^ = 4D/W§ = 252kHz > Fq 
as required. The Rayleigh range of the beam is ±6.3 cm 
from the cell center, so the beam expands by less than 8% 
within the cell, thus keeping a well-defined width along 
the cell. We adjusted the power of the light sheet to 150 
nW such that the power broadening is Tp = 20 Hz. 

In Fig. [H we plot the normalized measured dispersive 
spectrum of the Id beam configuration as a function of 
A, and the theory with the critical exponent as a single 
fit parameter. The data was logarithmically binned and 
averaged. No average was carried for detunings below 
300 Hz. Since we independently measured Fq it is not 
a fit parameter. The measured signal exhibits a narrow 
resonance feature with a width of Fq = 45 Hz, accom- 
panied with long tails ranging up to 8 kHz. We found 
P = 0.56 ± 0.01. We also reconstructed the absorptive 
part [m and observed clear deviations from a Lorentzian 
fit for the Id beam as shown in [9]. 

In Fig. O we present our main result. In this graph, 
we plot the absolute value of the complex EIT spec- 
trum derived from the measured dispersive part for the 
one-dimensional beam geometry (triangles) and for the 
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FIG. 3: Map of critical exponent (3 calculated from the exact 
theory [Eq. (5)] as a function of the diffusion broadening Td 
and the power broadening Fp in units of Fq. We extract P 
from the calculated spectrum using linear fit, with A in the 
range ±8 kHz. The linearity of the fit significantly increases 
towards zero and infinity. When Fp ^ oo at a fixed Td, P ^ 
0, signifying flat spectrum over the fit range. When the beam 
size tends to infinity, (Fd ^ 0) at a fixed Fp, the line-shape 
becomes a Lorentzian, and [3^1. The universal regime 
is where p 0.5. The triangle (circle) indicates the position 
where we performed our 1 — d (wide-area) measurements. The 
expected slope according to this prediction is 0.52, compared 
to the measured 0.56 d= 0.01. A similar calculation for the wide 
area beam yields 0.99, compared to the measured 0.97=b0.01. 



wide area beam (circles) as a function of the complex 
detuning size \s\ = y^A^ + (Fo/2)^ in units of Fo/2. In 
the universal regime, the complex spectrum is S{s) = 
s~^^P = 0.5. By writing s in polar representation, 
s = ^/S? + (ro/2)2e-^^^^*^"(2A/ro)^ ^-^^ dispersive part 

of the complex spectrum is sin[/3 arctan(2A/ro)]. In 
the experiment A was known and Fq was measured. If we 
let to be a single fit parameter and search for the best 
fit for the measured data (see Fig. 1), we obtain in the 
Id case jS = 0.56 ±0.01 over two decades in the frequency 
domain. For the wide area beam we obtain 0.97 ± 0.01 
over two decades. After obtaining an estimate for /3 we 
divide the measured spectrum by sin[/3 arctan(2A/Fo)] 
to obtain the absolute value of the complex spectrum 
[Fig. ([2])]. This should also be compared to the expected 
p = 0.52 and = 0.99 [Fig. (|3|)]. 

In Fig. [31 we plot the critical exponent as a function of 
the normalized diffusion broadening F d /^o and the nor- 
malized power broadening F^/Fq, which we calculated 
from the exact solution of Eq. (5). We verified that our 
chosen beam size and power are within the asymptotic 
regime where the universal theory is correct. Marked 
with a triangle is the expected critical exponent (0.52) 
for the parameters used in our experiment. For the 2d 
wide-area beam the expected critical exponent based on 



the exact solution [s*] is 0.99. In this case we verified that 
the presence of walls has a negligible effect on the critical 
exponent. We also studied the dependence on the opti- 
cal depth (OD). For OD < 10 the effect on the critical 
exponent is less than 0.5 percent. Higher values may al- 
ter the line-shape considerably due to the emergence of 
non-linear effects [l6| . 

The measured Laplace transform of the recurrence 
probability is in one-to-one correspondence with the FRT 
distribution in the time domain. In particular, when 
FRT(s) ~ 1 — 5^, and /3 < 1 as in our Id case, the 
asymptote of FRT(t) scales like t"^"^ 0. It follows that 
the mean first return time is divergent in the Id case. 
This stands in stark difference to the probability to stay 
in the dark presented in [17]. 

In conclusion, we found a regime where the EIT spec- 
trum has a universal form, which is the Laplace transform 
of the recurrence probability that depends only on the ef- 
fective dimensionality of the problem. Our result can be 
generalized to any kind of classical dynamics by replacing 
the diffusion propagator in Eq. (5) by the relevant center- 
of-mass propagator G(r, r', t) e.g. of anomalous diffusion, 
or of ballistic motion in billiards. In that case, following 
the same derivation it is found that the EIT spectrum in 
the universal regime is of the form G(0,0,5) = P(0,s). 
Staying in the realm of diffusion, our method can fa- 
cilitate measurements of interesting FRT and FPT dis- 
tributions in more complex settings such as diffusion in 
periodic and disordered (speckled) beam configurations 
and diffusion in the presence of traps. The latter are of 
particular interest, since the asymptotic survival prob- 
ability of the dark state in their presence is anomalous 
and has the form of a stretched exponential [5]. Simpler 
extensions are measuring the universal spectrum at 2d 
and measuring the Laplace transform of the first passage 
time from r = to r 7^ 0, which is proportional to the 
coherence at this point i?2i(r, 5). 
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